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ABSTRACT 

On 2008 January 10, the twin Solar Terrestrial Relations Observatory (STEREO) A and B 
spacecraft conducted a high time cadence study of the solar corona with the Extreme UltraViolet 
Imager (EUVI) instruments with the aim of investigating coronal dynamics. Observations of the 
three-dimensional propagation of waves within active region coronal loops and a measurement of 
the true coronal slow mode speed are obtained. Intensity oscillations with a period oi^l2 minutes 
are observed to propagate outwards from the base of a loop system, consistent with the slow 
magnetoacoustic mode. A novel analysis technique is applied to measure the wave phase velocity 
in the observations of the A and B spacecraft. These stereoscopic observations are used to infer 
the three-dimensional velocity vector of the wave propagation, with an inclination of 37 ± 6° to 
the local normal and a magnitude of 132 ± 9 and 132 ib 11 km s~^, giving the first measurement 
of the true coronal longitudinal slow mode speed, and an inferred temperature of 0.84 ± 12 MK 
and 0.84 ± 15 MK. 



Subject headings: MHD 
oscillations — Waves 

Introduction 
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The almost continual observations of the Solar 
and Heliospheric Observatory (SOHO) and Tran- 
sition Region and Coronal Explorer (TRACE) 
satellites have made observations of the solar 
corona readily available. During this time a mul- 
titude of coronal wave observations have been 
made, which may be interpreted as solutions to 
the Ma^neto-Hydrodynamic (MHD) wav e equa - 
tions as described by Edwi n fc RobertsI (Il983h . 
This recent consonance of theory and obser- 
vation has allowed the promise of using coro- 
nal seismology to derive hitherto indeterminable 
physical paramet ers of th e coro n al atmosphere 
as suggested by lUchidal (Il970h : [Roberts et al 
([l984). These principles have been used to 
estimate coronal pa rameters such as the mag- 
netic field strength ( Nakariakov fc OfmanI 2001 
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20021) . dissipative coe fficients (|Nakariakov et al.l 
1999!) , Alfven speeds (|Arregui et al 1 l2007h and 
density stratification ( Andries et al.ll2005[ ). 

A considerable limitation to current coronal 
seismology efforts is the unknown three dimen- 
sional geometry of coronal structures, since spa- 
tial measurements are confined to the two dimen- 
sional plane perpendicular to the observers line 
of sight. The STEREO mission consists of two 
spacecraft, one which orbits the Sun ahead of the 
Earth (STEREO A) and another that orbits be- 
hind (STEREO B). The STEREO spacecraft al- 
low the Sun to be observed from two different van- 
tage points, allowing stereoscopic measurements 
to be made, and the three dimensional nature of 
coronal structures to be inferred. 

Previous imaging observations of propagating 
slow-magnetoacoustic waves in coronal loops have 
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Fig. 1. — The 2048 x 512 pixel central strip of the Sun observed by EUVI during the high cadence observations, 
from the point of view of STEREO B (top) and STEREO A (bottom). The active region loops of interest are 
located at (390^', 10'') and (-350'', 10'') in heliocentric-cartesian coordinates relative to B and A respectively. 



been made using TRACE and SOHO/EIT. These 
are observed as intensity propagations in the coro- 
nal passbands with amplitudes on the order of a 
few percent and apparent propagation velo cities 
in the range 70-235 km s~^ (see Nightingale et al 



19991: iBerghmans fc Clett3ll999l : iRobbrecht et al 
2Q0ll : iDe Moort el et al.ll2002[ ). However, due to 
the inclination of the loop systems, these observed 
velocities are projections of the absolute speed per- 
pendicular to the line of sight. 

We present the first three dimensional observa- 
tions of coronal slow-magnetoacoustic wave prop- 
agation. These observations allow us to infer the 
three dimensional geometry of the wave propaga- 
tion vector and its magnitude, allowing the mea- 
surement of the true longitudinal slow mode speed 
and temperature within a coronal loop. 



2. Observations 

The observations were conducted on 2008 Jan- 
uary 10, as part of the Joint Observing Program 
(JOP) 200 - 'Multi-point, High Cadence EUV Ob- 
servations of the Dyna mic Solar Corona'. Using 
the EUVI instrument (jWuelser et al.ll2004l ). part 
of the Sun Earth Connection Cor onal and He- 
liospheric Investigation (SECCHI, iHoward et al 



[2008) instrument suite, to produce high cadence 
imaging of a sub-field region of interest. The 
aim of the STEREO/SECCHI led JOP was to 
use stereoscopic observations from the STEREO 
spacecraft, coupled to observations made in the 
Sun-Earth line, to produce a unique dataset exam- 
ining the three-dimensional, time-dependent, na- 
ture of the solar atmosphere. The observations 
consist of a series of 2048 x 512 pixel, sun-centered, 
sub-field images in the 171 A passband, with a 
fixed exposure and cadence of 4 s and 30 s, respec- 
tively. Table [1] lists relevant observational param- 
eters for both the STEREO A and B spacecraft. 
The high cadence observing program had a du- 
ration of two hours, beginning at 12:00:30UT for 
STEREO A and 12:00:50UT for STEREO B. The 
exposure timing is corrected for the light travel 
time between the two spacecraft, due to their dif- 
fering distance from the Sun, ensuring that obser- 
vations at both spacecraft locations originate from 
the Sun simultaneously. 

Within these data, we analyze the active region 
loop system whose footpoints are centered on solar 
x,y coordinates (-350", 10") in STEREO A and 
(390", 10") in STEREO B, where the coordinates 
are heliocentric-cartesian systems relative to the 
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Fig. 2.— The 200x200 pixel sub-field views of the coronal loop system from STEREO B (left) and STEREO 
A (right). The white lines outline the same coronal loop system, observed from differing points of view. 



STEREO A & B points of view. Figures □ and 
[2] show the relative location and geometry of the 
active region, as viewed from STEREO A & B, 
with a spacecraft separation of 44.5 degrees in the 
A-B-Sun plane. This loop system shows a complex 
pattern of outward intensity propagations, similar 
to that observed by previous EUV imagers; for the 
first time, however, we observe these propagations 
simultaneously from two different points of view. 
A novel analysis technique is applied to investigate 
the three dimensional nature of these propagations 
as outlined in the following section. 

3. Analysis 

3.1. Preparation of the data 

The EUVI images must be cleaned and cali- 
brated using a number of reduction procedures 
before the data can be analyzed. Using the SEC- 
CHLPREP routines available within the Solar- 
soft database, the data are corrected for detector 
bias, fiat field, and photometric calibration, which 
is applied to obtain the total number of counts 
per pixel. It is found that the series of images 
is affected by a systematic pointing 'jitter', pre- 
dominantly observed as a periodic variation of the 
pointing in the x direction. This pointing jitter is 
corrected by applying a 2D cross correlation tech- 
nique to long-lived, high contrast, features; in this 
case, stable plage structures observed at (335'V 



lO'O and (860'', 20'') in the A and B images respec- 
tively. The X pointing variation is found to have 
a sub-pixel amplitude on the order of a quarter of 
a pixel. 

A 200x200 pixel region of interest is selected, 
centered on the active region loop system in the 
A and B data. The solar rotation rate at the co- 
ordinates of the center of these sub-field regions is 
combined with the pointing jitter cross-correlation 
offsets and these combined offsets are used to co- 
align the images using cubic interpolation. Sub- 
sequently, two time series data cubes of 200x200 
pixels are then produced, centered on the active 
region loops, corrected and co-aligned as described 
above. 

3.2. Data analysis 

3.2.1. Time- distance analysis 

To investigate the oscillatory nature of the 
propagations traveling along the coronal loop sys- 
tem, the loop intensity is integrated along the path 
of the loops. The loops are outlined by a pair of 
arcs that are parallel to the visible loop strands 
and enclose the observed propagations along them. 
These arcs are then joined by successive cross- 
sections, two pixels wide, progressing along the 
length of the arcs and perpendicular to the loop 
strands. The intensity of the pixels along each 
cross-section is integrated and divided by the num- 
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Table 1 

Observational parameters of the EUVI high cadence observing program 



Observational Parameter 



stereo a 



stereo B 



Observations start 

FOV (pixels) 

Pixel scale (arcsec) 

Distance from Sun center (m) 

Image cadence (s) 

Exposure time (s) 

Passband (A) 

Heliographic Longitude (deg)^ 
Heliographic Latitude (deg)^ 



2008-01-10T12:00:30.004 
2048 X 512 
1.588 
1.4471540e+ll 
30.0 
4.0 
171 
21.295628 
-6.0627318 



2008-01-10T12:00:50.102 
2048 X 512 
1.590 
1.5073979e+ll 
30.0 
4.0 
171 
-22.987068 
-1.0544470 



^Heliographic longitude and latitude are Heliocentric Earth Equatorial (HEEQ), equiv- 
alent to the Stoneyhurst heliographic coordinates. 



Note. — During the observations the spacecraft had a separation of 44.5 degrees in the 
AB plane. 
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Fig. 3. — Time-distance images formed from the 
integrated intensity along the loops outlined in 
Fig El for STEREO B (left) and STEREO A 
(right). The abscissa indicates the observed dis- 
tance along the loop, perpendicular to the line of 
sight from each spacecraft to the photosphere. 



ber of cross-section pixels, to obtain the mean in- 
tegrated loop intensity as a function of distance 
along the loop system. Figure [2] shows the location 
of the defined arcs enclosing the loops for both the 
A and B images, where the central axis is parallel 
to the observed propagations. These integrated 
loop intensity profiles are then used to produce 
time-distance images of the intensity propagations 




Fig. 4. — Mean FFT of the integrated loop inten- 
sity from A (red) and B (blue), averaged over the 
lower ?^20,000 km as given in FigO 



along the magnetic field lines. 

Figure [3] shows the time-distance images for the 
coronal loops outlined in Fig [21 observed by the 
A and B spacecraft. Alternating light and dark 
bands of positive gradient indicate that a number 
of periodic outward intensity propagations are ob- 
served. They have a range of periodic timescales, 
most notably on the order of 700 s. Note the shal- 
lower gradient in STEREO A indicates a greater 
apparent velocity. 

Figure [4] indicates the average Fourier trans- 
form of the integrated loop intensity profiles for 
A and B over the lower ten cross-sections of the 
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Fig. 5. — Convolved time- distance images corre- 
sponding to Fig [3] filtered at 0.0014 Hz with a 
Gaussian filter of width 0.00014 Hz. 



time-distance images (or ?^20,000 km from the 
axes given in Fig [3]) . It is clear to see that the 
dominant oscillatory power is observed at a fre- 
quency of 0.0014 Hz from both spacecraft. There 
appears to be power around the 5-minute range 
(2.4-4 mHz) as m ight be expected for th is type 
of loop system (cf. IPe Moortel et al.ll2QQ2f ). How- 
ever, the observations from each spacecraft show 
differing power distributions. This suggests that 
the structures carrying these higher frequency os- 
cillations may be subject to line of sight effects 
from the point of view of each spacecraft, due to 
the differing temporal variability of the observed 
background plasma. 

3.2.2. Convolved time- distance analysis 

To investigate the spatial nature of the 0.0014 Hz 
oscillation in more detail, assuming stationary os- 
cillations, a convolution filtering technique is ap- 
plied to the time-distance images shown in Fig[3l 
A Gaussian filter of 0.14 mHz width is multiplied 
by the Fourier transform of the time series at 
each spatial location of the time-distance image. 
The position of the Gaussian filter is progressively 
scanned through the frequency range to isolate the 
oscillatory power at a particular frequency. The 
resulting signal to noise of the oscillation is in- 
creased allowing greater precision estimates of the 
propagation velocity to be made. Figure [5] shows 
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Fig. 6. — Least squares fit of the projected dis- 
tance along the loop system (perpendicular to the 
spacecraft line of sight) as a function of the prop- 
agation time from zero phase of the 0.0014 Hz os- 
cillation, for STEREO B (top) and STEREO A 
(bottom). 



the time-distance images with the Gaussian filter 
centered on 0.0014 Hz, thus isolating the spatial 
structure of the dominant propagating oscillation 
indicated in Fig[4l 

3.2.3. Propagation phase velocity 

To measure the apparent phase velocity of the 
0.0014 Hz propagating oscillation observed from 
STEREO A and B, the oscillation phase is de- 
termined as a function of distance using the con- 
volved time-distance images. The oscillation pa- 
rameters at each spatial position are calculated by 
applying the Bayesian spectral analysis code de- 
scribed in lMarsh et al.l ([2008). A single frequency 
oscillation model is applied to the data to deter- 
mine precise estimates of the oscillation parame- 
ters and their associated la uncertainties. The os- 
cillation phase propagation time is calculated from 
the oscillation phase, since the propagation time 
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from zero phase is given by r = ^/cj, where (j) and 
uj are the phase and angular frequency parameters 
determined by the Bayesian model. Regions where 
the phase jumps from to 27r are corrected by sub- 
tracting 27r, to provide continuity to the change in 
phase and phase propagation time. The phase ve- 
locity is then given by the inverse gradient from 
a least squares linear fit to the phase propagation 
time as a function of distance (Fig[6|). 

3.2.4- 3D wave propagation 

Once the propagating phase velocity is known 
from both the STEREO A and STEREO B points 
of view, then the three dimensional velocity vector 
may be calculated stereoscopically. The stereo- 
scopic 3D reconstruction software described in 
lAsch wanden et al.l (|2QQ8al lbl) is used to determine 
the three dimension al vector along; which the p rop- 
agation travels (see lAschwanden et aJ] l2QQ8bl . for 
a detailed description of the reduction process). 

Following the propagation direction, vectors 
within STEREO A and STEREO B are defined 
using four point splines, these are parallel to the 
coronal loops, and are along the central axis of 
the regions used to form the time-distance images 
(shown in the bottom panels of Fig.[7j). These 2D 
vectors in the STEREO A and B image planes are 
combined stereoscopically to constrain the unique 
solution of the true propagation vector in three di- 
mensions. The angle between the spacecraft line 
of sight and the 3D propagation vector, combined 
with the observed phase velocity, can be used to 
estimate the true propagation phase speed from 
each spacecraft. 

3.2.5. Error analysis 

The errors of the phase velocities measured in 
Fig. [6] are derived from the la errors on the filtered 
oscillation frequency and phase, returned by the 
Bayesian analysis code. These errors are propa- 
gated to determine the phase propagation time er- 
rors, which are applied in the least squares fitting 
routine used to determin e the propagation p hase 
velocity. As described in iMarsh et al.l (|2QQ8l ). the 



Bayesian code is able to return precise estimates 
of the oscillation parameters. It is the precision 
of the phase and frequency estimates, combined 
with the increased signal to noise of the convolved 
time-distance analysis, that allows precise mea- 



surements of the propagation phase velocity. 

The Aschwanden stereoscopy code determines 
an error on the three dimensional loop coordinates 
as the error along the line of sight in the reference 
frame of STEREO A. lAschwanden et al.l |2008b) 
estimates these height errors assuming that the 
position of a loop in the EUVI A & B images can 
be selected to within one pixel by an automated 
loop fitting procedure. Due to the extended struc- 
ture of the observed loop system, we estimate the 
loop vectors defined within A and B have a two 
pixel uncertainty. The error in the height of the 
points that define the three dimensional propaga- 
tion vector are used to estimate the error range 
of the propagation inclination. The dependence 
on cosine results in an asymmetric error in the 
inclination. This error in the propagation inclina- 
tion is propagated with the error in the measured 
phase velocity from A and B, to determine the er- 
ror in the true propagation phase speed for each 
spacecraft. 

4. Results 

4.1. 0.0014 Hz wave propagation 

Figure [3] shows the time-distance images formed 
from the loops outlined in Fig. O The complex 
pattern of periodic propagations indicates that a 
number of frequencies propagate upwards along 
the magnetic field of the loops. Figure (H shows 
that the propagations are dominated by a fre- 
quency at 0.0014 Hz. The convolved time-distance 
images in Fig. [5] present the time-distance data iso- 
lated for the 0.0014 Hz frequency. It can been 
seen from Fig. [5] that the 0.0014 Hz oscillation 
propagates with a constant velocity along the loop 
when observed from both spacecraft, but that 
the propagation in the STEREO A observations 
has a shallower gradient than that observed from 
STEREO B. This suggests that the propagation 
has a greater apparent velocity perpendicular to 
the line of sight when viewed from STEREO A 
compared to STEREO B. Considering the obser- 
vations of the loop system in Fig. [21 this may 
be expected as the loop system appears to be 
more inclined away from the observers line of sight 
when viewed from STEREO A, thus increasing the 
velocity component perpendicular to the line of 
sight. 



6 




Fig. 7. — Derived 3D geometry of the wave propagation along the coronal loops. Bottom panels show the 
observed region from STEREO A & B overplotted with the projected propagation vector (solid line), and 
the arcs used to form the running difference analysis from Fig. [2] (dotted line). The top panels show the 
projected geometry of the propagation vector with the view rotated 90 degrees about the x-axis looking 
towards solar north. 



4.2. Observed phase velocity 

Figure [6] shows the projected distance along 
the coronal loop system against the phase prop- 
agation time of the 0.0014 Hz oscillation. The 
gradient gives the phase propagation velocity in 
the plane perpendicular to the observation line 
of sight, where the negative value defines upward 
propagation from the base of the loop system. The 
data is fitted with a linear function (there is no 
significant acceleration trend as a function of dis- 
tance along the loop), thus implying a constant 
phase velocity. As the time-distance images sug- 
gest, the phase velocity measured from STEREO 
A (104.1 ±0.2 km s~^) is greater than that mea- 
sured from STEREO B (62.9 ± 0.1 km s"^), in- 
dicating that the propagation has a greater in- 
clination to the line of sight when observed from 
STEREO A. 

4.3. 3D wave propagation 

Using: the stereoscopic reconstruction code 
( Aschwanden et al ] [2008b), the suggested prop- 
agation inclination relative to the A and B space- 



craft is verified, with the propagation inclination 
of 52.lt^;f relative to the STEREO A line of 
sight, an inclination of 28.4^2^° relative to the 
line of sight of STEREO B, and an inclination 
of 36.8^5"]^ to the local solar normal. The true 
propagation speed derived from the phase velocity 
observed from STEREO A is 132.0t^;^ km s-^ 
Using the phase velocity observed from STEREO 
B, the derived true propagation speed is 1?>2.2'^^^ 
km s~^. Table [2] lists the observed velocities, rel- 
ative geometry, and the absolute phase speed of 
the propagating wave. 

5. Discussion 

The results presented here represent the first 
measurements of the true propagation speed of the 
longitudinal slow mode within a coronal loop. The 
propagation speed of the slow mode can be derived 
from the dispersion relation f or waves trapped 
within a magnetic cylinder ( Edwin fc RobertsI 
19831 ). The slow mode propagates within a cylin- 



drical geometry with a phase speed q < Vph < Cs • 
In the long wavelength limit, in comparison with 
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Table 2 

Parameters of the propagating wave 



Propagation Parameter STEREO A STEREO B 

Measured phase velocity (km s~^) 104.1 ± 0.2 62.9 ±0.1 

Inclination to LOS (degrees) 52.ltl\t 28Atl\t 

Inclination to local normal (degrees) 36.8'^^'\ 36.8'^^'\ 

Reconstructed phase speed (km s~^) 132. olg ^ 132.2^1%^ 



the diameter of the structure, the slow mode prop- 
agates at the tube speed c^, and in the short wave- 
length limit propagates at the sound speed Cs . The 
waves observed here tend to the long wavelength 
limit and, strictly, propagate at the tube speed 



Cf = 



(1) 



where Ua is the Alfven speed. In the case where 
Va > Cs, then Ct = Cg. 

If we assume a low-beta coronal plasma with 
a typical Alfven speed of 1000 km s~^ and the 
measured tube speed of 132 km s~^, the tube and 
sound speed differ by less than 1%; thus we may 
approximate the tube speed to be equal to the 
sound speed. The sound speed for a weakly bound 
collisional plasma such as the low corona may be 
approximated using the ideal gas equation of state, 
to derive the temperature. 



(2) 




T = 



cliirrip 



(3) 



where the adiabatic index 7 = 5/3, p is the plasma 
pressure, p is the plasma density, k is the Boltz- 
mann constant, T is the temperature, nip is the 
proton mass, and assuming a fully ionized hydro- 
gen and helium plasma of coronal abundances, the 
mean molecular weight /i = 1.27 and the mean 
particle number Nj^ = 1.92. 

The measured phase speeds of 132.0lg"5 km s~^ 
and 132.2^^^4^ km s-\ from STEREO A and B 
respectively, are used to derive infer the tempera- 
ture of OMtoAi MK and 0.84t^;J? MK. This de- 
rived temperature is very close to the peak of the 
EUVI 171 A response functions, as shown in Fig.[8l 
The 171 A bandpass is dominated by emission from 
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Fig. 8. — Response functions of the EUVI 171A 
channel for STEREO A (solid) and STEREO B 
(dashed). The vertical lines indicate the loop tem- 
perature (solid) and errors (dotted), inferred from 
the STEREO A phase speed results. 

the Fe IX 171 A line formed around these temper- 
atures, suggesting that the waves are observed by 
emission from Fe IX 171 A. The correspondence 
between the sound speed at the peak response 
temperature and the measured phase speed gives 
further evidence that the slow magnetoacoustic 
mode is observed. 

It should be remembered, however, that these 
results are based on the assumption of low plasma 
beta. It is conceivable that over-dense structures 
may exist in the corona, in which case the Alfven 
speed can be of the same order of magnitude as 
the sound speed. The tube and sound speeds may 
then be significantly different, from Equation [TJ 
With accurate spectroscopic estimates of the tem- 
perature and density, this allows the possibility of 
discriminating between the tube and sound speeds 
for seismological purposes, such as determining 
the magnetic field strength, or adiabatic index, for 
example; this will be addressed in a following pa- 
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per. 

The large variation in the observed phase ve- 
locities from STEREO A & B demonstrates the 
significance of the projection effect on previous 
observational studies of longitudinal modes us- 
ingEIT and TRACE, such as De Moortel et all 
( 2QQ2[ ). The projected phase velocity in the plane 
of the sky cannot be approximated to the phase 
speed and used to derive meaningful physical pa- 
rameters. Using single point observations, it is 
not possible to disentangle the projection effect to 
study the physical causes of any variation in ob- 
served phase speeds. 

The results for the phase speed derived from 
both the STEREO A and STEREO B observa- 
tions agree very well within the errors. This close 
agreement between the results suggests that the 
correct 3D propagation geometry has been deter- 
mined and validates the implementation of the 
Aschwanden et ahl ( 2Q08bl ) stereoscopic code. 



The constant velocity shown in Figs. [5] and [6l 
combined with the fact that the propagation fol- 
lows a linear path from the stereoscopic results, 
suggests that a constant adiabatic index and equi- 
librium temperature is maintained along the loop 
strands, probably due to the high thermal conduc- 
tion along the magnetic field. 

There has been a limited number of previ- 
ous observations of low frequency oscillations 
of t his order i n cor onal loops. For example, 
[rbrra das e t all (|2QQ4 ) find an :^0.0016 Hz oscil- 
lation in TRACE o bservations of fan like loops, 
[Mcintosh et_alJ ([2008) also observe a 0.0015 Hz 
oscillation in TRACE loops associated with plage. 
It is possible that lar g er sca le studies such as that 
of iDe M oortel et all (}2Q02') could fail to detect 
these low frequency oscillations due to a selection 
effect of using a wavelet analysis and time series 
of limited duration, since wavelet analysis has a 
detection cutoff at low frequencies, due to the 
cone of influence of the wavelet transform. Also, 
the use of running difference analysis over a time 
range designed to detect periodicities in the range 
3-5 minutes will selectively exclude other periods 
and detection of lower frequencies. 

6. Conclusions 

The measurement of the phase speed of the lon- 
gitudinal slow magnetoacoustic mode in coronal 



loops is presented for the first time. A 0.0014 Hz 
intensity oscillation propagating within a coronal 
loop system is observed by both EUVI instruments 
on the STEREO A and B spacecraft. High pre- 
cision estimates of the wave phase velocity, ob- 
served in the plane of the sky from each spacecraft, 
are determined using convolved time-distance im- 
ages and a Bayesian oscillation code. The absolute 
phase speed is then calculated using the three di- 
mensional geometry of the propagation obtained 
from the stereoscopic observations. 

The measured phase speeds of 132. Olg ^ km s~^ 
and 132.2^3^4^ km s~^ are consistent with the 
expected speed of the slow magnetoacoustic 
mode. The inferred loop strand temperature of 
OMtoll MK and OMtofi MK, under the as- 
sumption of low plasma beta, is close to the peak 
of the 171A response function. Thus, the observed 
phase speed corresponds to the sound speed at 
this temperature, further strengthening the slow 
magnetoacoustic mode interpretation. 

In principle, it is possible to use the propagating 
slow mode to measure the coronal magnetic field 
strength. Temperature and density diagnostic ob- 
servations are necessary to determine if the slow 
mode tube speed can be discriminated from the 
sound speed within the errors. This possibility is 
currently under investigation and will be reported 
in a following paper. 
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